Effects of Bending on Raman-active Vibration Modes of Carbon Nanotubes 
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We investigate vibration modes and their Raman activity of single-walled carbon nanotubes that 
are bent within their intrinsic elastic limits. By implementing novel boundary conditions for density- 
functional based tight-binding, and using non-resonant bond polarization theory, we discover that 
Raman activity can be induced by bending. Depending on the degree of bending, high-energy 
Raman peaks change their positions and intensities significantly. These effects can be explained by 
migration of nodes and antinodes along tube circumference. We discuss the challenge of associating 
the predicted spectral changes with experimental observations. 
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Technological applications of carbon nanotubes 
(CNTs) are based on their exceptional mechanical and 
electronic properties. Due to advances in fabrication 
and manipulation, CNTs have become one of the most 
prominent building blocks for nanoscale materials design. 
Their electronic and mechanical properties can be used 
in numerous applications for nanoclcctronics, hydrogen 
and energy storage material, sensors, or high-strength 
composites i 

Among mechanical properties 2 -, vibrations arc rel- 
evant in heat dissipation 3 -, sensors^ and nanotube 
identification 5 -. Due to the large number of applications, 
vibrations have been investigated extensively. In particu- 
lar, most experimental studies use Raman spectroscopy^, 
a method that is able to achieve even single nanotube 
resolution^ 

In practice, because CNTs are long, they bend. Bend- 
ing is observed in isolated CNTs between electrodes^, 
or in "paper"—, "forests"—, ringsii, and composite 
systems^ made out of CNTs. It appears that bending 
is ubiquitous in experiments — and challenging to study 
theoretically. Most previous theoretical Raman studies 
are for straight tubes, because modeling of bent systems 
has been computationally too expensive^. Modeling of 
bending with classical methods is straightforward, but 
has been used to study force moments and strains^ 3 -, 
bucklin g 14 : 15 and other large-scale mechanical properties 
that result from rather high curvature 3 

In this work we investigate how vibrations and Ra- 
man spectra are affected when CNTs are bent slightly, 
within their intrinsic elastic limits. To accomplish this, 
novel boundary conditions are introduced which allow 
quantum mechanical modeling of bending with compu- 
tationally feasible system size. We show that vibrations 
undergo systematic changes that significantly alter the 
high-frequency Raman spectra. The spectral changes can 
be understood via simple physical principles. 

We use density-functional based tight-binding (DFTB) 
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FIG. 1: Schematic of the boundary condition. Wave func- 
tions satisfy ^(rg) = V(^-±0( r s)), where yb is a point on 
either boundary (with one k-point along tube axis); opera- 
tion 7Z±^(rj) rotates Yj an angle ±<j> around the wedge apex. 
We define the vector Vij from atom i to atom j as equal to 
shortest of vectors (TZ^(rj) — rj), (lZ-#(rj) — r 4 ) or (rj — n), 
where r$ and Vj are atom positions within the simulation box. 
For some atom pairs crosses the boundary and renders 
Newton's third law invalid. 



mcthoc— i to calculate forces, optimize systems 1 ^, and 
calculate vibrational cigenmodes. The method has been 
used successfully for vibrational analysis of carbon nan- 
otubes, also related to Raman activity! 18 : 19 Raman spec- 
tra are calculated by non-resonant bond polarization 
metho d 20 ' 21 . This method has some restrictions and lim- 
its direct comparison between resonant experimcnt a 22 ' 23 , 
but suffices for the scope of this paper. The details of 
our approach are given in Ref. l24l . 

To model bent nanotubes quantum mechanically, we 
introduce novel "periodic wedge boundary conditions" 
where the CNT appears as a slice of a torus, as shown in 
Fig. [TJ Some complications arise from fixed quantization 
axis and require mild approximations 3?. Henceforth all 
directions refer to the fixed Cartesian coordinates shown 
in Fig. [TJ 

We stress that in our calculations the bending of CNTs 
is only due to boundary conditions. All atoms are free to 
move in the unit cell; no constraints are applied (those 
would cause severe artifacts for the vibration modes). 
Note that the length of the tube is automatically op- 




FIG. 2: Raman spectra of bent (13,0), (6, 6) and (13,4) carbon nanotubes in different polarization pictures. Polarization picture 
is defined by directions of incident and scattered light as shown by symbols on the right. Bending increases linearly from zero 
(upmost lines) to O = 4.2 % for (13,0), = 2.6 % for (6,6), and = 4.8 % for (13,4) tubes (lowest lines; corresponds to wedge 
angle <j) = 10° for all tubes). The symbols with arrows refer to symmetry of modes for straight tubes in line group notation^. 
Peaks were Lorenzian broadened with full width at half maximum of 5 cm -1 . Small symbols in the spectra of (6,6) tube with 
zz-polarization refer to the subscripts of original mode symmetries. 



timizcd since the atoms can freely move in the radial 
direction of the unit cell. 

The single- walled CNTs we investigate are: semi- 
conducting (13,0) tube with 10.2 A diameter and 
21.3 A length, metallic (6,6) tube with 8.1 A diame- 
ter and 27.1 A length, and metallic (13,4) tube with 
12.1 A diameter and 21.9 A length. To have a common 
measure for the degree of bending for tubes with differ- 
ent chiralities and diameters, we define a dimcnsionlcss 
variable 

6 = D/2R, 

where D is tube's diameter and R the radius of cur- 
vature measured from tube axis. Hence G = 1.0 (100 
%) corresponds to maximum bending (torus with a van- 
ishing hole); buckling of a nanotube takes place above 
O ~ 10 % 14 i 15 and the range for bondings in experiments 
is estimated 9 = 0.05 % ... 5 which is the range 

under our focus. 

Our main results are embedded in Fig. [21 showing 
Raman spectra of high-energy modes by systematically 
varying O. Low-energy modes are not shown because 
they are insensitive to bending, with respect to both en- 
ergy and intensity. For example, the radial breathing 



mode is left nearly intact because bending mostly affects 
bonds parallel to tube axis. On the contrary, within G- 
band -related high-energy modes bending induces sys- 
tematic changes: (i) emergence of new peaks, (ii) inten- 
sity reductions of "original" (straight-tube) peaks, (Hi) 
significant energy shifts and (iv) splitting of peaks into 
smaller ones. In (13,0) tube E+ mode even deactivates 
for moderate bending. The rich spectra of chiral (13,4) 
tube is further enriched while bent, but still shows similar 
systematics to achiral tubes. 

We begin analyzing these results by looking at what 
happens in the spatial structure of the vibration modes 
as tubes are bent. To visualize this, we show a qualita- 
tive view of the nodes and antinodes for selected modes 
in Fig. [31 Two-dimensional E modes within the same po- 
larization picture are similar and uniform in the direction 
of the tube, differing only in the number of nodes along 
tube circumference. In xy and zz polarizations the am- 
plitudes are in (13,0) tube along tube axis and in (6,6) 
tube along tube circumference (for xz vice versa). As 
the tube is bent, the nodal structure starts migrating to- 
wards the outer or inner side of the torus. Amplitudes are 
also affected, but the modes remain uniform along tube 
axis and the transition is smooth. In other words, modes 
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FIG. 3: (color online) Qualitative view of the nodes and antin- 
odes of selected vibration modes for (13,0) and (6,6) tubes. 
Thick blue line represents an antinode, red an antinode with 
opposite phase and vanishing line a node; the actual direction 
of the vibration amplitude is either circumferential or along 
tube axis. Modes are for straight or bent tubes as indicated 
on the left. The symbols below show the polarization pictures 
where the bent modes become most visible. We were unable 
to identify one mode for (6,6) tube (faint symbol). 



in slightly bent tubes can be always identified with the 
symmetric modes in straight tubes. 

A qualitative change for bent tubes, as shown in Fig. [4] 
as a function of 0, is the change in bond lengths and bond 
stiffnessess. The bonds in inner parts shorten, bonds 
in outer parts lengthen, with up to 5 % variation for 
= 4%. Because carbon bonds in honeycomb structure 
are stiff, the result shows a surprisingly strong (~ 50 %) 
variation in the effective (nearest neighbor) spring con- 
stants between the inner and outer parts of the tube. 

These observations help us to understand the energy 
shifts, remembering that lj ~ \fk. Let us take the zz- 
polarization spectrum of (6,6) tube in Fig. [2] as an ex- 
ample. The energies of original Aq and E% symmetry 
modes increase; it is because the antinodes concentrate 
in the inner part of the torus where the spring constants 
are larger (see Fig. [3]). The energies of vibrations orig- 
inating from Ei and E$ symmetry modes decrease; it 
is because the antinodes concentrate in the outer part 
of the torus where the spring constants arc smaller (sec 
Fig. [3}. Furthermore, looking at Fig. [3] one can realize 
that the energy of E^ mode decreases more than E$ be- 
cause the antinodes of E± concentrate in outer part more 
strongly. This is part of a more general trend observ- 
able in Fig. [3] modes with many circumferential nodes 
are less affected because antinodes are more equally dis- 
tributed as tubes are bent. The same observations apply 
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FIG. 4: Variation of bond lengths for (13,0) (left) and (6,6) 
(right) tubes as a function of bending parameter Q. Bonds are 
divided into mostly parallel (||) and mostly perpendicular (_L) 
bonds with respect to tube axis (in and out refer to the sides 
of torus). Inset shows the effective nearest neighbor spring 
constant, calculated by scaling single graphene layer; the in- 
set was used to map the bond length scale (left-hand axis) 
into spring constant scale (right-hand axis). The equilibrium 
spring constant 970 N/m and thermal expansion coefficient 
5.13 • 10~ 6 K _1 , used to make the inset, are in overall agree- 
ment with reference results . 27 ' 28 



to other spectra. Note that in Fig. [3] bending modifies 
two-dimensional modes pairwise the same way (antin- 
odes migrate in same direction), and bending does not 
lift the degeneracy of the modes because energy is for 
both modes either increased or decreased. Hence the ap- 
pearance of peak splittings in Fig. [2] is due to originally 
different modes, not due to lifted degeneracy. 

Why does bending cause migration of nodes and antin- 
odes along the circumference? Consider a simple model: 
calculate the eigenmodes of a simple, one-dimensional 
ring of atoms connected by harmonic springs, where 
spring constants are modulated such that on one side 
of the ring they are larger, on the other side smaller. For 
this modified system the symmetric eigenmodes change 
so that nodes and antinodes concentrate either on the 
region of strong bonds, or on the region of weak bonds, 
depending on the mode in question. In nanotubes the 
modes and amplitudes are three-dimensional, but the ba- 
sic mechanism remains the same. 

Finally, let us investigate the induced Raman activity 
of Fig. [21 Consider for example the E^ mode of (6,6) 
tube in Fig. [3] that is originally Raman inactive with xy 
and zz polarizations. Bending breaks the symmetry in 
x-direction and causes E^ mode to resemble A^ and E^ 
modes — modes that arc Raman active for straight tubes. 
Due to this resemblance Ef and Aq get nearly equal Ra- 
man intensities^. More generally, from y-antisymmetric 
modes, that can be excited by y-polarized light but that 
have broken x-symmetry, some will "leak" the polariza- 
tion into cc-direction, making the mode active with xy 
polarization (e.g. left-hand E^'s for (6,6) tube in Fig. [3]). 
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Further, from y-symmetric modes, that have broken x- 
symmetry, some become active in zz polarization (e.g. 
right-hand E^'s for (6,6) tube in Fig. [3]). For (13,4) 
tube similar principles apply, but corresponding analysis 
is somewhat more complicated due to the chirality of the 
tube. 

The arguments above are biased towards the non- 
resonant bond polarizability model, but we stress that 
most arguments are related to vibrational eigcnmodes 
(peak positions, node migration), and are independent of 
the method to calculate the Raman intensities. There- 
fore most observations should be consistent with more 
complete theories. 

In fact, even the most complete theory would have 
problems with experimental interpretation. Spectra are 
calculated only for a piece of potentially complex curved 
CNT system. We confirmed that the tiny x-component, 
that tube axis has near boundaries, is not the origin for 
Raman activity "leakage" between polarization pictures. 
But if nanotube slice is a part of more complete torus, 
situation becomes more complicated and direct compar- 
ison less sensible, because tube axis mixes with other 
directions. Because the situation depends crucially on 
the experimental setups, we cannot make here general 
interpretations. 

There have been few experiments aiming for direct ob- 
servation of Raman spectra for bent nanotube a 11 ' 29 . In 
Ref. [29| the G band peak was observed to broaden and 
shift lower in energy, which was attributed to the in- 
creased bond lengths in bent tubes. Because for ideally 



bent tubes some modes should also decrease in energy, it 
is likely that in this experiment the tubes were not only 
bent but also stretched; shift is due to stretching and 
broadening due to bending. Our highest bending limit, 
= 5%, was derived from microscopic images^ 2 -, but 
high density of such bendings must involve defects. It is 
because a bent tube must be anchored via mechanical or 
chemical bonds; a rough estimate for the bending energy 
of a tube with any chirality is -Ehend ~ 20 -6 2 (eV/atom), 
which for 9 = 5% bending requires significant 5 eV of 
anchoring energy per 100 atoms in the tube. Therefore 
defect-free extreme bending should occur only in singu- 
lar parts of compound structures, and the Raman inten- 
sity from these parts is expected to remain comparably 
small. On the other hand, Raman measurements for sin- 
gle CNTs could show these effects visibly, especially un- 
der experimentally feasible controlled bending^. 

To conclude, this work provides understanding into ef- 
fects that bending causes for vibrational spectra of CNTs. 
The effects are significant, but can be systematically ex- 
plained with general principles, which should be valid 
for multi-walled nanotubes or even for non-carbon nan- 
otubes. 
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